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Based on 5.8 x 1Q 7 J/tp events detected in BESII, the branching fractions of J/tp — > <j>rj and (prj 
are measured for different r\ and rj' decay modes. The results are significantly higher than previous 
measurements. An upper limit on B(J/tp —* <frw°) is also obtained. 



PACS numbers: 13.25. Gv, 12.38.Qk, 14.40.Gx 

I. INTRODUCTION 

The decay of the J / tp into a vector and pseudoscalar 
meson pair, J/tp — ► VP with V and P represent- 
ing vector and pseudoscalar mesons, can proceed via 
strong and electromagnetic reactions. A well mea- 
sured set of all possible decays of J/tp — > VP al- 
lows one to systematically study the quark gluon con- 
tents of pseudoscalar mesons, SU(3) breaking, as well 
as determine the electromagnetic and doubly sup- 
pressed OZI amplitudes in two-body J/tp decays [1]. 
MARKIII [2, 3] and DM2 [4] measured many J/tp -> 
VP decays and obtained the 77 — 77' mixing angle, the 
quark content of the 77 and 77', and much more. 

Recently, a sample of 5.8 x 10 7 J/ip events was 
accumulated with the upgraded Beijing Spectrome- 
ter (BESII) [5], which offers a unique opportunity 
to measure precisely the full set of J/tp — > VP de- 
cays. In an earlier analysis based on this data set, 
the branching fraction of J/ip — > 7r + 7r~7r° was mea- 
sured to be (2.10 ± 0.12)% [6], which is higher than 
the PDG [7] value by about 30%. This indicates a 
higher branching fraction for J/tp — > pir than those 
from older experiments [8], since the dominant dy- 
namics in J/tp — > 7T + 7r~7r° is J/tp — > pir. Therefore, 
remeasuring the branching fractions of all J/tp — > VP 
decay modes becomes very important. In this paper, 
J/tp — > 07T° , 077, and (pij' are studied, based on the 
BESII 5.8 x 10 7 J/tP events. 

II. THE BES DETECTOR 

The upgraded Beijing Spectrometer detector (BE- 
SII) is located at the Beijing Electron-Positron Col- 
lider (BEPC). BESII is a large solid-angle mag- 
netic spectrometer which is described in detail in 
Ref. [5]. The momentum of charged particles is de- 
termined by a 40-layer cylindrical main drift cham- 
ber (MDC) which has a momentum resolution of 
CT p /p=1.78% v /l + p 2 (p in GeV/c). Particle iden- 



tification is accomplished using specific ionization 
(dE/dx) measurements in the drift chamber and time- 
of-flight (TOF) information in a barrel-like array of 
48 scintillation counters. The dE/dx resolution is 
° 'dE/dx — 8.0%; the TOF resolution for Bhabha events 
is <jtof — 180 ps. Radially outside of the time-of- 
flight counters is a 12-radiation- length barrel shower 
counter (BSC) comprised of gas proportional tubes 
interleaved with lead sheets. The BSC measures the 
energy and direction of photons with resolutions of 
a E /E ~ 21%-/E (E in GeV), cr = 7.9 mrad, and 
g z = 2.3 cm. The iron flux return of the magnet is 
instrumented with three double layers of proportional 
counters (MUC) that are used to identify muons. 

A GEANT3 based Monte Carlo package (SIMBES) 
with detailed consideration of the detector perfor- 
mance is used. The consistency between data and 
Monte Carlo has been carefully checked in many high 
purity physics channels, and the agreement is reason- 
able. The detection efficiency and mass resolution for 
each decay mode are obtained from a Monte Carlo 
simulation which takes into account the angular dis- 
tributions appropriate for the different final states [9] . 

III. ANALYSIS 

In this analysis, the <p meson is observed in its 
K + K~ decay mode, and the pseudoscalar mesons 
are detected in the modes: tt° — > 77; 77 — > 77, 
77r + 7T~, and it + it~it ; and 77' — > 77, 77r + 7r~(7p), 
and 7r + 7r~77 (77 — ► 77). Using multiple 77 and 77' 
decay modes allows us to crosscheck our measure- 
ments, as well as obtain higher precision. Possible 
final states of J/tp — > (pn , <prj, and <pr/' are then 
K + K~jj, K + K~tt + it~j, and K + K~ir + ir~jj. Can- 
didate events are required to satisfy the following com- 
mon selection criteria: 

1. The events must have the correct number of 
charged tracks with net charge zero. Each track 
must be well fitted to a helix, originating from 



3 



the interaction region of R xy <0.02 m and \z\ < 
0.2 m, and have a polar angle, 8, in the range 
| cos 6» | < 0.8. 

2. Events should have at least the minimum num- 
ber of isolated photons associated with the dif- 
ferent final states. Isolated photons are those 
that have energy deposited in the BSC greater 
than 60 MeV, the angle between the direction at 
the first hit layer of the BSC and the developing 
direction of the cluster less than 30°, and the 
angle between photons and any charged tracks 
larger than 10°. 

3. For each charged track in an event, XpidW i s 
determined using both dE/dx and TOF infor- 
mation: 



A. J/xp — > </>77 

Events with two oppositely charged tracks and at 
least two or three isolated photons are selected. A 
4C-fit is performed to the K + K~jj hypothesis, and 
X 2 < 15 is required. To reject possible background 
events from J/tp — > jK + K~tt° , the 4C-fit probability 
for the assignment J/tp — > K + K~-f-f must be larger 
than that of K + K~jjj. 

After this selection, the scatter plot (Figure 1) of 
ttik+k- versus m 77 shows two clusters corresponding 
to <jyq' and (prj, but there is no clear accumulation of 
events for <fiir° . To obtain the distribution recoil- 
ing against </>, the K + K~ invariant mass is required 
to be in the (p mass region, \m K+K - — 1.02| < 0.02 
GeV/c 2 . 



X 2 pid (i)=X 2 dE /d x M+Xtof (*) 



A charged track is identified as a ir or K if its 
Xpid 1S l ess than those for any other assignment. 
To reject background events, two charged tracks 
are required to be identified as kaons in J/tp — > 
<pn°. For the other channels, at least one charged 
track must be identified as a kaon in the event 
selection. 

4. The selected events arc fitted kincmatically. The 
kinematic fit adjusts the track energy and mo- 
mentum within the measured errors so as to sat- 
isfy energy and momentum conservation for the 
given event hypothesis. This improves resolu- 
tion, selects the correct charged-particle assign- 
ment for the tracks, and reduces background. 
When the number of photons in an event ex- 
ceeds the minimum, all combinations are tried, 
and the combination with the smallest \ 2 is re- 
tained. 

The branching fraction is calculated using 



B{J/tp -» <j>P) 



N, 



obs 



Nj/t,, -e-B(cP^ K+K-) ■ B(P -> X) ' 

where N bs is the number of events observed (or the 
upper limit), Nj/^ is the number of J/ip events, 
(5.77 ± 0.27) x 10 7 , determined from the number of 
inclusive 4-prong hadronic decays [10], e is the detec- 
tion efficiency obtained from Monte Carlo simulation, 
and B((p -> K+K~) and B(P -» X) arc the branch- 
ing fractions of — > K + K~ and pseudoscalar decays 
from the PDG [7], respectively. 




m(YY) (GeV/c 2 ) 



FIG. 1: Scatter plot of m K + K - versus m 77 for J/ip 
K + K~^ events. 



1. J/tp-* (pTT 

Figure 2(a) shows the m 77 invariant mass distribu- 
tion after the above selection; no clear tt° signal is ob- 
served. The Bayesian method is used to determine the 
upper limit on the J/tp — > (pir° branching fraction. A 
Brcit-Wigner convoluted with a Gaussian plus a poly- 
nomial background function are used to fit the m 77 
spectrum. The ir° mass and width are fixed to PDG 
values. The mass resolution, obtained from Monte 
Carlo simulation, is 17.7 MeV/c 2 . At the 90% confi- 
dence level, the number of (pir° events is 24. Taking 
into account the detection efficiency (16.63 ± 0.20)%, 
the upper limit on the branching fraction is 

B(J/i/> -» </>7T°) < 5.10 x 10~ 6 
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3. J/ip -> 0r/ 
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FIG. 2: The invariant mass distribution of m 77 for J/ip — > 
077 events. The curves are the results of the fit described 
in the text. 



2. J/ip -» 0?7 

Figure 2(b) shows the distribution; an rj signal 
is clearly seen. The fit of this distribution with a Breit- 
Wigner convoluted with a Gaussian plus a second or- 
der polynomial background function gives 2086 ± 42 
cp-q events with a r\ mass of 549.0 ± 0.5 MeV/c 2 . The 
background events, 152 ± 17, are estimated from the cp 
sidebands, defined by 0.98 GeV/c 2 < m K+K - < 1.00 
GeV/c 2 and 1.04 GeV/c 2 < m K+K - < 1.06 GeV/c 2 . 
After subtracting background and correcting for de- 
tection efficiency, (19.98 ± 0.22)%, the J/ip -> cp-q 
branching fraction is obtained 

B(J/tp -> cjyq) = (8.67 ± 0.19) x 10~ 4 , 

where the error is statistical only. 



The distribution of m 77 in rf mass region recoil- 
ing against the cp is shown in Figure 2(c). A fit of 
the rj' peak with a Breit-Wigner and a second order 
backgound polynomial yields 68 ± 15 cprj' events with 
the peak at 958.1 ± 2.6 MeV/c 2 . No obvious signal is 
observed for the distribution of recoiling against 
sidebands (0.98 GeV/c 2 < m K+K - < 1.0 GeV/c 2 
and 1.04 GeV/c 2 < m K+K - < 1.06 GeV/c 2 ). The 
detection efficiency is (18.57 ± 0.22)%, and the corre- 
sponding branching fraction is determined to be 

B{J/tP -> cprj') = (6.10 ± 1.34) x 10~ 4 , 

where the error is only the statistical error. 



B. J/ip — > 4> r yK + -K 

For J/ip — > cjyq, r\ — > 77r + 7r~, events with four well- 
reconstructed charged tracks and at least one isolated 
photon are required. To select the pions and kaons 
from amongst the tracks, 4C fits are applied for one 
of the following three cases: (1) if only one charged 
track is identified as a kaon using particle identifica- 
tion, then the other charged tracks are assumed, one 
at a time, to be a kaon, while the other two are as- 
sumed to be pions; (2) if two charged tracks are iden- 
tified as kaons, then the other two tracks are assumed 
to be pions; (3) if three or four charged tracks are 
identified as kaons, then the particle identification in- 
formation is ignored and all combinations of two kaon 
and two pion tracks are kinematicaly fitted. For each 
case, the hypothesis with the smallest x 2 is selected. 
We further require that the probability of the 4C fit for 
the J/ip — > K + K~Tt + Tt~"/ assignment is larger than 
those of K + K~ir + ir~ and A" + A"~7r + 7r~77. 

is shown 



The scatter plot of m K + K - versus m 



in Figure 3, where J/ip — > <pr] and J/ip — > cjyq' decays 
are apparent. For the scatter plot of m w + T - versus 
m 7T + T -, shown in Figure 4, the tf - p signal cor- 



responds to the decay rj' —> 
is from r\ — > 77r + 7r~ and r\ 
events. 



7/9. The other cluster 
background 



7T 7T + 7r 



1. J/ip — > 4>q 

Figure 5 shows the 77r + 7r~ invariant mass recoil- 
ing against the 0, defined by \m K + K - — 1 .02 1 < 0.02 
GeV/c 2 . A clear -q signal is observed. The peak on 
the left side of the 77 in Figure 5 comes from J/ip — > 
cpr] (rj — ► 7T + 7r~7r°) with one photon missing; this is 
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m(K + K) (GeV/cf) 



FIG. 3: Scatter plot of m K + K - versus m 77r+7r - for J/ip — ► 
K + K~ tt + it~ 7 events. The band below the rj signal comes 
from rj — + 7r + 7r~7r events. 



branching fraction is 

B(J/r/> -» = (9.79 ± 1.02) x 10~ 4 , 
where the error is only the statistical error. 
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FIG. 5: Distribution of m 7 ^.+ w 



7T+7T 7 



for J/i/> - 

events. Dots with error bars are data, and the curves are 
the results of the fit described in the text. 
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FIG. 4: Scatter plot of m„+ n - versus m 77r + n - for J/^i — ► 
4>"yn + -K~ events. The rj' - p signal corresponds to the decay 
rf — > 7p. The other cluster is from rj — > 77r + 7r~ and 
r/ — > 7r°7r + 7r _ background events. 



confirmed by Monte-Carlo simulation. This peak can- 
not be described by a simple Brcit-Wigncr due to its 
asymmetric shape. To obtain the shape of the peak, 
a Monte-carlo sample of J/ip — > <^>?7 (?7 — > 7r + 7r~7r°) is 
generated and a fit is made to the peak. The jit + tt~ 
mass distribution is then fitted with this shape, a 
Brcit-Wigner to describe the rj signal, and a polyno- 
mial background. The fit, shown in Figure 5, yields 
134 ± 14 rj events with a mass at 548.9 ± 0.9 MeV/c 2 . 
The detection efficiency obtained from Monte Carlo 
simulation is (10.32 ± 0.16)%, and the corresponding 



2. J/tp -> <j)7j' 



After requiring \m K+K - - 1.02| < 0.02 GeV/c 2 
and 0.3 GeV/c 2 < m n+n - < 0.95 GeV/c 2 , the dis- 
tribution of 77r + 7r~ invariant mass recoiling against 
the <j> is shown in Figure 6; a fit with a Breit- 
Wigner convoluted with a Gaussian and a second or- 
der polynomial gives 462 ± 29 events with a peak at 
957.4±0.7 MeV/c 2 . The detection efficiency obtained 
from Monte Carlo simulation is (9.80 ± 0.16)%, and 
the branching fraction obtained is 

B(J/4> -» <p7j') = (5.64 ± 0.35) x 10" 4 . 



C. J ftp — » (f)-K + ir 77 

For the rj — > 7r + 7r~7r° case, events with four well 
reconstructed charged tracks and at least two iso- 
lated photons are selected. A 4C kinematic fit to the 
_fT + ii~7r + 7r~77 hypothsis is applied, as described in 
Section IIIB for J/ip — > 077T+7T - , and the case with 
the smallest x 2 is selected. 

After the above selection and with the requirement 
that m 77 be consistent with a 7r°, (0.095 GeV/c 2 
< m 77 < 0.175 GeV/c 2 ), the J/ip — > (f>tj decay is 
clearly observed in the scatter plot of m K + K - ver- 
sus m Tr + Tr - 77 , shown in Figure 7(a). Requiring 0.5 
GeV/c 2 < m 77 < 0.6 GeV/c 2 , the scatter plot in 
Figure 7(b) shows clean <frq' signals. The decays of 
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m(yJiV) (GeV/c 



FIG. 6: The distribution of m ylv + n - for events of the type 
J/tp — > (j>py; the curves are the result of the fit described 
in the text. 



rj — ► 7r + 7r _ 7T° and 77' — > 7r + 7r~?7 are also observed in 
the scatter plot of m 77 versus m„.+ „.-,,, shown in Fig- 
ure 8. 



The m K + K - invariant mass spectrum recoiling 
against the 77, shown in Figure 9, is used to get the 
(j>rj signals. A Breit-Wigner convoluted with a Gaus- 
sian to account for the <j> mass resolution plus a sec- 
ond order polynomial are used to fit the m K + K - 
mass distribution. A total of 350±11 events with a 
cj> mass at 1020.4 ± 0.3 MeV/c 2 from (p decay are ob- 
tained in the fit, which using the detection efficiency 
of (5.81 ± 0.12)% corresponds to a branching fraction 
of 

B{ J/ip -> ^77) = (9.41 ± 0.30) x 10~ 4 . 
Here, the error is only the statistical error. 



2. J/tp -> (pr, 1 

After requiring 0.5 < to 77 < 0.6 GeV/c 2 and 
m T +j- < 0.45 GeV/c 2 , the 7r + 7r~77 mass recoiling 
against the 4>{\m K + K - — 1.02| < 0.02 GeV/c 2 ), shows 
a clean 77' peak, as seen in Figure 10. No clear signal is 
observed for <f> sidebands (0.98 GeV/c 2 < m K + K - < 
1.0 GeV/c 2 and 1.04 GeV/c 2 < m K+K - < 1.06 
GeV/c 2 ). The fit of m T+7r - 77 yields 198 ± 12 events 
with a peak at 959.2 ± 1.4 MeV/c 2 , and the detection 
efficiency for this channel is (7.83±0.14)%, which gives 




m(K + K") (GeV/cT) 



FIG. 7: Scatter plots for m K + K - versus m n + n -^o and 
m K+K - versus m 7T + 7T - v for J/tp — > K + K~ tt + -k~ 77 
events. 
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FIG. 8: Scatter plot of m 77 versus m 7r+7r - 77 for J/tp 
07r + 7r~77 candidate events. 
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1. MDC 



efficiency 
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1.15 



m(K + K") (GeV/c 2 ) 



FIG. 9: The m K + K - distribution for J/ip — > 
K + K _ 7r + 7r _ 7r° events. The curves are the results of the 
fit described in the text. 



B(J/ip -> <jrq') = (5.11 ± 0.31) x 10~ 4 . 
Here, the error is statistical only. 




m(7iVr)) (GeV/c ) 



FIG. 10: The m^+ 7r - 77 distribution for J/tp — > <f>TT + n^n 
candidate events. The curves are the result of the fit de- 
scribed in the text. 



D. Systematic Errors 

In this analysis, the systematic errors on the 
branching fractions mainly come from the following 
sources: 



The MDC tracking efficiency is measured in 
clean channels like J/ip — > AA and ip(2S) — > 
7r + 7r~ J/tp, J/tp — > /i + 'fi~ . It is found that the Monte 
Carlo simulation agrees with data within 1-2% for 
each charged track. Therefore 4% is taken as the sys- 
tematic error on the tracking efficiency for the chan- 
nels with two charged tracks and 8% for the channels 
with four charged tracks in the final states. 



2. Particle ID 

The particle identification (PID) efficiency of the 
kaon is studied from J/tp — > K + K~ir° and J/tp — + <pi]. 
The results indicate that the kaon PID efficiency for 
data agrees well with that of the Monte Carlo sim- 
ulation in the kaon momentum region less than 1.0 
GeV/c. In the analysis of J/tp — > <pir°, where two 
charged tracks are required to be kaons, the PID effi- 
ciency difference between data and Monte Carlo sim- 
ulation is about 3.4%. In other decay modes, at least 
one charged track is required to be identified as a kaon, 
so the difference from PID is less than 1%. Here, the 
difference of the PID efficiencies between data and 
Monte Carlo simulation is taken as one of the system- 
atic errors. 



3. Photon detection efficiency 

For the decay modes analyzed in this paper, one or 
two photons are involved in the final states. The pho- 
ton detection efficiency is studied from J /if) — ► p°n° 
in Ref. [6]. The results indicate that the difference 
between the detection efficiency of data and MC sim- 
ulation is less than 2% for each photon. 



4- Kinematic fit 

The kinematic fit is a useful tool to improve res- 
olution and reduce background. The systematic er- 
ror from the kinematic fit is studied with the clean 
channel J/ip — > tt + tt~tt , as described in Ref. [6]. 
The conclusion is that the kinematic fit efficiency dif- 
ference between data and Monte Carlo simulation is 
about 4.1%. Using the same method, the decay mode 
J/ip — * 7r + 7r~7T + 7r~7r° is also analyzed, and the kine- 
matic fit efficiency difference between data and Monte- 
Carlo is about 4.3%. In this paper, 5% is conserva- 
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tivcly taken to be the systematic error from the kine- 
matic fit for all analyzed decay modes. 

5. Selection criteria 

The systematic errors for additional selection cri- 
teria in specific decay modes are estimated by com- 
paring the efficiency difference with and without the 
criterion or replacing it with a very loose requirement. 
The study indicates that they are not large compared 
with other systematic errors. The results are listed in 
Table 1 

6. Uncertainty from hadronic interaction model 

Different simulations of the hadronic interaction 
lead to different efficiencies. In this analysis, two mod- 
els, FLUKA [11] and GCALOR [12], are used in simu- 
lating hadronic interactions in the Monte-Carlo. The 
difference of the detection efficiencies from these two 
Monte Carlo models is about 3%, which is taken as 
the systematic error. 

7. Uncertainty of background 

The uncertainties of the background in each channel 
are estimated by changing the background shape in 
the fit. The results are listed in Table I. 

8. Intermediate decay branching fractions 

The branching fractions of <p — > K + K~ and the 
pseudoscalar decays arc taken from the PDG. The er- 
rors of these branching fractions are systematic errors 
in our measurements and are listed in Table I. 

The systematic error contributions studied above, 
the error due to the uncertainty of the number of J/ip 
events, and the statistical error of the Monte-Carlo 
samples are all listed in Table I. The total systematic 
error is the sum of them added in quadrature. 



IV. RESULTS AND DISCUSSION 

The branching fractions of J/ip decaying into 4>tt°, 
<prj, and (fyq', measured into different final states, are 
listed in Table II. The average value is the weighted 
mean of the results from the different decay modes, 
and the PDG value is the world average taken from 
Ref. [7]. The world averages mainly come from 
Marklll and DM2. The results obtained here are not 
in good agreement with previous measurements. Just 
as for the branching fraction of J/ip — > 7r + 7r~7r , the 
branching fraction of J/ip — > (p>r\ and <jnf are higher 
than those in the PDG. 

In this paper, we measured the branching fractions 
of J/ip decays into <p plus a pseudoscalar. The three 
branching fractions arc not sufficient for a detailed 
study of pseudoscalar mixing, SU(3) breaking, and the 
contribution from doubly suppressed OZI processes 
using the phenomenological model in Ref. [1]. How- 
ever the inconsistency between the results from BE- 
SII and those from former measurements emphasize 
the importance for such a study. After measuring the 
other decay modes of J/ip VP, such as J/ip — > ivn , 
U7], tor]', prj, pi]', and K*K, it will be important to 
extract physics with all the relevant measurements 
again. 
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